A genome-wide association study provides insights into the genetic etiology of 57 essential and non-essential trace elements in humans.

Trace elements are important for human health but may exert toxic or adverse effects. Mechanisms of uptake, distribution, metabolism, and excretion are partly under genetic control but have not yet been extensively mapped. Here we report a comprehensive multi-element genome-wide association study of 57 essential and non-essential trace elements. We perform genome-wide association meta-analyses of 14 trace elements in up to 6564 Scandinavian whole blood samples, and genome-wide association studies of 43 trace elements in up to 2819 samples measured only in the Trøndelag Health Study (HUNT). We identify 11 novel genetic loci associated with blood concentrations of arsenic, cadmium, manganese, selenium, and zinc in genome-wide association meta-analyses. In HUNT, several genome-wide significant loci are also indicated for other trace elements. Using two-sample Mendelian randomization, we find several indications of weak to moderate effects on health outcomes, the most precise being a weak harmful effect of increased zinc on prostate cancer. However, independent validation is needed. Our current understanding of trace element-associated genetic variants may help establish consequences of trace elements on human health.

Large scale plasma proteomics identifies novel proteins and protein networks associated with heart failure development.

Heart failure (HF) causes substantial morbidity and mortality but its pathobiology is incompletely understood. The proteome is a promising intermediate phenotype for discovery of novel mechanisms. We measured 4877 plasma proteins in 13,900 HF-free individuals across three analysis sets with diverse age, geography, and HF ascertainment to identify circulating proteins and protein networks associated with HF development. Parallel analyses in Atherosclerosis Risk in Communities study participants in mid-life and late-life and in Trøndelag Health Study participants identified 37 proteins consistently associated with incident HF independent of traditional risk factors. Mendelian randomization supported causal effects of 10 on HF, HF risk factors, or left ventricular size and function, including matricellular (e.g. SPON1, MFAP4), senescence-associated (FSTL3, IGFBP7), and inflammatory (SVEP1, CCL15, ITIH3) proteins. Protein co-regulation network analyses identified 5 modules associated with HF risk, two of which were influenced by genetic variants that implicated trans hotspots within the VTN and CFH genes.

Insight into the relationship between resting heart rate and atrial fibrillation: a Mendelian randomization study.

AIMS: A low resting heart rate (RHR) implies a more efficient heart function and a lower risk of cardiovascular disease. However, observational studies have reported a U-shaped association between RHR and atrial fibrillation (AF). In contrast, Mendelian randomization (MR) studies have found an inverse causal association between RHR and AF. Hence, the causal nature of the relationship is not clear. The aim is to investigate the causal association and its shape between RHR on AF using linear and non-linear MR (NLMR). METHODS AND RESULTS: Linear and non-linear MR were performed on individual-level data in the Trøndelag Health Study (HUNT) and UK Biobank (UKB). HUNT consists of 69 155 individuals with 7,062 AF cases, while UKB provides data on 431 852 individuals with 20 452 AF cases. The linear MR found an inverse relationship between RHR and AF with an OR = 0.95 [95% confidence interval (CI): 0.93-0.98] and OR = 0.96 (95% CI: 0.95-0.97) per unit decrease in RHR in HUNT and UKB, respectively. The NLMR was supportive of an inverse linear relationship in both HUNT and UKB for RHR values <90 beats per minute (bpm). Several sensitivity analyses were also consistent. CONCLUSION: In contrast with the current observational knowledge of RHR and AF, an inverse causal association between RHR and AF was demonstrated in both linear and non-linear MR for RHR values up to 90 bpm. Further exploring the underlying mechanisms of the genetic instrument for RHR may shed light on whether pleiotropy is biasing this association.

Educational attainment, health outcomes and mortality: a within-sibship Mendelian randomization study.

BACKGROUND: Previous Mendelian randomization (MR) studies using population samples (population MR) have provided evidence for beneficial effects of educational attainment on health outcomes in adulthood. However, estimates from these studies may have been susceptible to bias from population stratification, assortative mating and indirect genetic effects due to unadjusted parental genotypes. MR using genetic association estimates derived from within-sibship models (within-sibship MR) can avoid these potential biases because genetic differences between siblings are due to random segregation at meiosis. METHODS: Applying both population and within-sibship MR, we estimated the effects of genetic liability to educational attainment on body mass index (BMI), cigarette smoking, systolic blood pressure (SBP) and all-cause mortality. MR analyses used individual-level data on 72 932 siblings from UK Biobank and the Norwegian HUNT study, and summary-level data from a within-sibship Genome-wide Association Study including >140 000 individuals. RESULTS: Both population and within-sibship MR estimates provided evidence that educational attainment decreased BMI, cigarette smoking and SBP. Genetic variant-outcome associations attenuated in the within-sibship model, but genetic variant-educational attainment associations also attenuated to a similar extent. Thus, within-sibship and population MR estimates were largely consistent. The within-sibship MR estimate of education on mortality was imprecise but consistent with a putative effect. CONCLUSIONS: These results provide evidence of beneficial individual-level effects of education (or liability to education) on adulthood health, independently of potential demographic and family-level confounders.

Lowering of Circulating Sclerostin May Increase Risk of Atherosclerosis and Its Risk Factors: Evidence From a Genome-Wide Association Meta-Analysis Followed by Mendelian Randomization.

OBJECTIVE: In this study, we aimed to establish the causal effects of lowering sclerostin, target of the antiosteoporosis drug romosozumab, on atherosclerosis and its risk factors. METHODS: A genome-wide association study meta-analysis was performed of circulating sclerostin levels in 33,961 European individuals. Mendelian randomization (MR) was used to predict the causal effects of sclerostin lowering on 15 atherosclerosis-related diseases and risk factors. RESULTS: We found that 18 conditionally independent variants were associated with circulating sclerostin. Of these, 1 cis signal in SOST and 3 trans signals in B4GALNT3, RIN3, and SERPINA1 regions showed directionally opposite signals for sclerostin levels and estimated bone mineral density. Variants with these 4 regions were selected as genetic instruments. MR using 5 correlated cis-SNPs suggested that lower sclerostin increased the risk of type 2 diabetes mellitus (DM) (odds ratio [OR] 1.32 [95% confidence interval (95% CI) 1.03-1.69]) and myocardial infarction (MI) (OR 1.35 [95% CI 1.01-1.79]); sclerostin lowering was also suggested to increase the extent of coronary artery calcification (CAC) (ß = 0.24 [95% CI 0.02-0.45]). MR using both cis and trans instruments suggested that lower sclerostin increased hypertension risk (OR 1.09 [95% CI 1.04-1.15]), but otherwise had attenuated effects. CONCLUSION: This study provides genetic evidence to suggest that lower levels of sclerostin may increase the risk of hypertension, type 2 DM, MI, and the extent of CAC. Taken together, these findings underscore the requirement for strategies to mitigate potential adverse effects of romosozumab treatment on atherosclerosis and its related risk factors.

Dyslipidemia and Risk of Preeclampsia: A Multiancestry Mendelian Randomization Study.

BACKGROUND: Preeclampsia is a leading cause of maternal morbidity, and dyslipidemia has been associated with preeclampsia in observational studies. We use Mendelian randomization analyses to estimate the association between lipid levels, their pharmacological targets, and the risk of preeclampsia in 4 ancestry groups. METHODS: We extracted uncorrelated (R2<0.001) single-nucleotide polymorphisms strongly associated (P<5×10-8) with LDL-C (low-density lipoprotein cholesterol), HDL-C (high-density lipoprotein cholesterol), and triglycerides from genome-wide association studies of European, admixed African, Latino, and East Asian ancestry participants. Genetic associations with risk of preeclampsia were extracted from studies of the same ancestry groups. Inverse-variance weighted analyses were performed separately for each ancestry group before they were meta-analyzed. Sensitivity analyses were conducted to evaluate bias due to genetic pleiotropy, demography, and indirect genetic effects. RESULTS: The meta-analysis across 4 ancestry groups included 1.5 million subjects with lipid measurements, 7425 subjects with preeclampsia, and 239 290 without preeclampsia. Increasing HDL-C was associated with reduced risk of preeclampsia (odds ratio, 0.84 [95% CI, 0.74-0.94]; P=0.004; per SD increase in HDL-C), which was consistent across sensitivity analyses. We also observed cholesteryl ester transfer protein inhibition-a drug target that increases HDL-C-may have a protective effect. We observed no consistent effect of LDL-C or triglycerides on the risk of preeclampsia. CONCLUSIONS: We observed a protective effect of elevated HDL-C on risk of preeclampsia. Our findings align with the lack of effect in trials of LDL-C modifying drugs but suggest that HDL-C may be a new target for screening and intervention.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis.

BACKGROUND: Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. RESULTS: To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. CONCLUSIONS: Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals.

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects.

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Taller height and risk of coronary heart disease and cancer: A within-sibship Mendelian randomization study.

Background: Taller people have a lower risk of coronary heart disease but a higher risk of many cancers. Mendelian randomization (MR) studies in unrelated individuals (population MR) have suggested that these relationships are potentially causal. However, population MR studies are sensitive to demography (population stratification, assortative mating) and familial (indirect genetic) effects. Methods: In this study, we performed within-sibship MR analyses using 78,988 siblings, a design robust against demography and indirect genetic effects of parents. For comparison, we also applied population MR and estimated associations with measured height. Results: Within-sibship MR estimated that 1 SD taller height lowers the odds of coronary heart disease by 14% (95% CI: 3-23%) but increases the odds of cancer by 18% (95% CI: 3-34%), highly consistent with population MR and height-disease association estimates. There was some evidence that taller height reduces systolic blood pressure and low-density lipoprotein cholesterol, which may mediate some of the protective effects of taller height on coronary heart disease risk. Conclusions: For the first time, we have demonstrated that the purported effects of height on adulthood disease risk are unlikely to be explained by demographic or familial factors, and so likely reflect an individual-level causal effect. Disentangling the mechanisms via which height affects disease risk may improve the understanding of the etiologies of atherosclerosis and carcinogenesis. Funding: This project was conducted by researchers at the MRC Integrative Epidemiology Unit (MC_UU_00011/1) and also supported by a Norwegian Research Council Grant number 295989.

Trans-ethnic Mendelian-randomization study reveals causal relationships between cardiometabolic factors and chronic kidney disease.

BACKGROUND: This study was to systematically test whether previously reported risk factors for chronic kidney disease (CKD) are causally related to CKD in European and East Asian ancestries using Mendelian randomization. METHODS: A total of 45 risk factors with genetic data in European ancestry and 17 risk factors in East Asian participants were identified as exposures from PubMed. We defined the CKD by clinical diagnosis or by estimated glomerular filtration rate of <60 ml/min/1.73 m2. Ultimately, 51 672 CKD cases and 958 102 controls of European ancestry from CKDGen, UK Biobank and HUNT, and 13 093 CKD cases and 238 118 controls of East Asian ancestry from Biobank Japan, China Kadoorie Biobank and Japan-Kidney-Biobank/ToMMo were included. RESULTS: Eight risk factors showed reliable evidence of causal effects on CKD in Europeans, including genetically predicted body mass index (BMI), hypertension, systolic blood pressure, high-density lipoprotein cholesterol, apolipoprotein A-I, lipoprotein(a), type 2 diabetes (T2D) and nephrolithiasis. In East Asians, BMI, T2D and nephrolithiasis showed evidence of causality on CKD. In two independent replication analyses, we observed that increased hypertension risk showed reliable evidence of a causal effect on increasing CKD risk in Europeans but in contrast showed a null effect in East Asians. Although liability to T2D showed consistent effects on CKD, the effects of glycaemic phenotypes on CKD were weak. Non-linear Mendelian randomization indicated a threshold relationship between genetically predicted BMI and CKD, with increased risk at BMI of >25 kg/m2. CONCLUSIONS: Eight cardiometabolic risk factors showed causal effects on CKD in Europeans and three of them showed causality in East Asians, providing insights into the design of future interventions to reduce the burden of CKD.

Global Biobank Meta-analysis Initiative: Powering genetic discovery across human disease.

Biobanks facilitate genome-wide association studies (GWASs), which have mapped genomic loci across a range of human diseases and traits. However, most biobanks are primarily composed of individuals of European ancestry. We introduce the Global Biobank Meta-analysis Initiative (GBMI)-a collaborative network of 23 biobanks from 4 continents representing more than 2.2 million consented individuals with genetic data linked to electronic health records. GBMI meta-analyzes summary statistics from GWASs generated using harmonized genotypes and phenotypes from member biobanks for 14 exemplar diseases and endpoints. This strategy validates that GWASs conducted in diverse biobanks can be integrated despite heterogeneity in case definitions, recruitment strategies, and baseline characteristics. This collaborative effort improves GWAS power for diseases, benefits understudied diseases, and improves risk prediction while also enabling the nomination of disease genes and drug candidates by incorporating gene and protein expression data and providing insight into the underlying biology of human diseases and traits.

The HUNT study: A population-based cohort for genetic research.

The Trøndelag Health Study (HUNT) is a population-based cohort of ∼229,000 individuals recruited in four waves beginning in 1984 in Trøndelag County, Norway. Approximately 88,000 of these individuals have available genetic data from array genotyping. HUNT participants were recruited during four community-based recruitment waves and provided information on health-related behaviors, self-reported diagnoses, family history of disease, and underwent physical examinations. Linkage via the Norwegian personal identification number integrates digitized health care information from doctor visits and national health registries including death, cancer and prescription registries. Genome-wide association studies of HUNT participants have provided insights into the mechanism of cardiovascular, metabolic, osteoporotic, and liver-related diseases, among others. Unique features of this cohort that facilitate research include nearly 40 years of longitudinal follow-up in a motivated and well-educated population, family data, comprehensive phenotyping, and broad availability of DNA, RNA, urine, fecal, plasma, and serum samples.

Risk of lower respiratory tract infections: a genome-wide association study with Mendelian randomization analysis in three independent European populations.

OBJECTIVE: Lower respiratory tract infections (LRTIs) are a leading cause of morbidity and mortality worldwide. Few studies have previously investigated genetic susceptibility and potential risk factors for LRTI. METHODS: We used data from the UK Biobank, Trøndelag Health Study (HUNT), and FinnGen to conduct a genome-wide association study (GWAS). Cases were subjects hospitalized with LRTI, and controls were subjects with no such hospitalization. We conducted stratification and interaction analyses to evaluate whether the genetic effect of LRTI differed by sex or smoking. Mendelian randomization (MR) analyses were conducted to identify the unconfounded relationship between cardiometabolic risk factors and LRTI. RESULTS: A total of 25 320 cases and 575 294 controls were included. The 15q25.1 locus reached genome-wide significance in the meta-analysis (rs10519203: OR 0.94, p 3.87e-11). The protective effect of effect allele of rs10519203 was present among smokers (OR 0.90, 95%CI 0.87-0.92, p 1.38e-15) but not among never-smokers (OR 1.01, 95%CI 0.97-1.06, p 5.20e-01). In MR analyses, we found that increasing body mass index (OR 1.31, 95%CI 1.24-1.40, p 3.78e-18), lifetime smoking (OR 2.83, 95%CI 2.34-3.42, p 6.56e-27), and systolic blood pressure robustly increased the risk of LRTIs (OR 1.11, 95%CI 1.02-1.22, p 1.48e-02). CONCLUSION: A region in 15q25.1 was strongly associated with LRTI susceptibility. Reduction in the prevalence of smoking, overweight, obesity, and hypertension may reduce the disease burden of LRTIs.

The power of genetic diversity in genome-wide association studies of lipids.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.